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ABSTRACT 

The quark-hadron phase transition in the early universe can produce inho- 
mogeneities in the distribution of nucleons, which in turn affect the primordial 
nucleosynthesis. In all the investigations of this problem it has been assumed that 
the degree of supercooling of the quark-gluon plasma after the phase transition 
is large enough to produce a significant rate of nucleation of hadrons. Using the 
latest results of finite temperature lattice QCD and the finite size scaling theory, 
we argue that the degree of supercooling is in fact extremely small and hence the 
nucleation rate is negligible. 
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The implications of the transition from the deconfined quark-gluon plasma 
phase to the confined hadron phase in the early universe for the primordial nucle- 
osynthesis have been studied extensively in recent years [1-11]- The following phase 
separation scenario has been used in these investigations. In the early universe 
when the temperature is high the unconfined quarks and gluons form a plasma. 
As the universe cools due to expansion and reaches the phase transition tempera- 
ture, Tc, the confined hadron phase does not appear immediately since the phase 
transition is assumed to be of first order. The universe has to cool somewhat 
below Tc before the probability for nucleating bubbles becomes significant. The 
supercooled phase is metastable and decays into hadronic bubbles. 

The rate of nucleation of hadron bubbles has been calculated [4,12] using the 
theory of homogeneous nucleation [13]. In this theory the new phase is assumed to 
arise through spontaneous fluctuations in the supercooled phase. The nucleation 
rate per unit volume, p(T), at a temperature, T, is given by 

p{r^)^CT^eM-^^j^) (1) 

where C is a constant of order unity, u, the surface tension of the hadronic bubble, 
L, the latent heat per unit volume of the transition and 77, the degree of supercool- 
ing, ?7 = (Tc — T)/Tc. As the supercooling is assumed to be small, after a hadron 
bubble has been nucleated it grows explosively as a deflagration bubble [12,14] 
and drives a weak shock-wave which expands into the surrounding quark phase 
with a velocity just above sound speed Vs — l/-\/3 . The plasma gets reheated by 
the shock- wave and further nucleation is inhibited due to the steep fall-off in Eq. 
(1). The fraction of the universe affected by the shock front preceding the growing 
bubble is 

g{t) = j\{t')^vl{t-t'fdt' , (2) 

where tc is the time when the universe cools through Tg. Here the temperature T 
in p{T) has been replaced by time t using the Einstein relation between T and the 
age of the universe,t, 

*~ 27r(4l7r)i/2 

In calculating t the quark degrees of freedom have been included^. The total 



^ We have used here the ideal gas expression for the energy density. While it 
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number density of nucleation sites is then 

Nn= / f{t)p{t)dt , (4) 

J to 

where 

f{t) = l-g{t) . (5) 

Since the supercoohng, that is, 77 is small, p{t) rises rapidly with t and therefore 
f{t) behaves almost like a step function d{th — t) [4, 12], where is the time when 
the universe has been reheated. A^^ can then be approximated as 

Nn= r P{t)dt . (6) 

Jto 

By expanding p{t) around th one obtains [4] , 

pit) ^ p{th) ew{-Oi{th - t)) , (7) 

where 

2 / rp \ 2 

a = 1.6528 X 10-^^1 — —] , (8) 

7^ is a dimensionless combination of cr, L, and Tg, 
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- c 

and rjh can be calculated from the condition f{th) — 0, i.e.. 



dtp{t')-v',{t-t')^l . (10) 



Substituting Eq. (7) in Eq. (10) and carrying out the integration one finds the 
equation determining ijfi to be, 

a;^ exp(a;) = 2.21 x 10"^ (^ ^^^^ ^ , (11) 

where 

16n 7^ 



X = 



3 Vt 



is now known[15] that near Tg there are strong non-perturbative effects, they will 
affect Eq. (3) by a factor of 0(1) and have therefore been ignored. 



(We have put C = 1.). An approximate solution of Eq. (11) given in ref. 10 
is r]h = 0.47. Note that the RHS of Eq. (11) is totally dominated by the very 
large ratio of the Planck scale and the QCD scale, which is typically ~ 10^^. 
Thus, a large variation of 7, over several orders of magnitude but still sufficiently 
small compared to the ratio, changes the solution very little, making it almost 
independent of the details of the QCD thermodynamics. 

One of the parameters important for primordial nucleosynthesis is the mean 
separation per nucleation site, I ~ Nn , which can be obtained from Eq. (6) by 
evaluating the integral using Eqs. (7)-(8) and the solution to Eq. (11) given above. 
One obtains 

;=(87r)V3!!£ . 

a 

A more detailed calculation by Meyer et al [10] gives 

1 = 4.38— . (12) 

a 

During the subsequent evolution of the universe, when the hadronization of 
the quark-gluon plasma is complete, the mean free paths of neutrons and pro- 
tons in the surrounding medium of hadrons, electrons, photons and neutrinos 
become different leading to a segregation such that high-density proton-rich and 
low-density neutron-rich regions develop. Nucleosynthesis then takes place among 
these inhomogeneously distributed nucleons [3]. It is clear that for this picture of 
nucleosynthesis to be valid / should be greater than the comoving diffusion length 
of the protons, which is ~ 0.5 m. The fluctuations would otherwise damp out 
entirely before the onset of nucleosynthesis. 

It is apparent from Eq. (8), (9) and (12) that the degree of supercooling and 
hence / depend only on the bulk thermodynamic properties, namely, Vs, cr, L and 
Tc of the strongly interacting matter. One knows, however, from finite size scaling 
theory [16] that the degree of metastability is also decided by the size of the system 
relative to its correlation length, ^(Tc) at the transition point. Indeed , the width, 
AT of the metastable region for a first order phase transition goes to zero as the 
volume, V, increases: 

^.^-1 = ^, (13) 

where A is a constant of 0(1) if V is measured in units of ^{Tc). 

In order to proceed further and see the consequences of Eq. (13) we shall 
assume that the world of quenched quantum chromodynamics (QCD) is a good 
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approximation to ours. In particular, we wish to use the precise information avail- 
able from lattice QCD in this approximation. Furthermore, we shall assume that 
asymptotic scaling is valid so that we can relate the lattice results to the continuum 
theory. We shall return to the discussion of the effects of these assumptions later. 
The best estimates from lattice QCD for the quantities required are as follows [15, 
17, 18]: 

L/T^ = 1.8 ±0.1 

a/T^ = 0.024 ± 0.004 . (14) 
Tc = 235 ± 42 MeV 

Note that the transition temperature is much larger here than the value used in 
Ref. 4 from the MIT bag model. The other parameters are roughly the same. 
In fact, using Tg as a scale, one sees that it is really the latent heat,L, which is 
much smaller here ( by almost a factor of 8), while the surface tension is within 
the range used in Ref. 4. Using these values in Eqs.(9) and (12), one obtains 
7 = 2.066 X 10~^ and I = 2.84 ± 1.25 cm, which is much less than what is required 
even if we ignore the question of the size of the universe at T = Tc. 

If we now ask how much supercooling we can expect by taking size of the 
universe into account the situation becomes considerably worse. We know from 
lattice QCD [19] that i) the physical correlation length in the deconfined phase 
is ^{Tc)Tc = 1.38 ± 0.24 and ii) the width of the metastable region in the QCD- 
coupling space for a system of volume V is given by 

-^3s 0.984±0.039 

v) ' . (15) 



A/? = 0.29 



Further, by assuming the validity of the asymptotic scaling relation, the width in 
the coupling /3 can be related to the width of the metastable region in temperature: 
AT/Tc — 1.122A/3. Using the expression for the age of the universe, Eq. (3), the 
volume of the universe at temperature is found to be 

^r|^7.1><10-(H?i^)\ (16) 

Substituting in Eq. (15), one obtains the degree of supercooling possible in the 
early universe to be 

A.rp / T' \ 2.952±0.117 

— =r7 = 9.28xl0-^6±2.2/_lc \ . 

Tr ' \200MeV ) ^ ' 
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Thus for any value of in the range given in Eq. (14) the universe practically 
never supercools at all. Even the miniscule amount indicated by Eq. (17) will not 
be sufficient for any bubble nucleation, since the nucleation probability goes to 
zero as exp(— 1/?7^) as 77 — > 0. 

Let us now discuss whether our assumptions made earlier could have led to 
the negative result. Two crucial inputs were the validity of the asymptotic scaling 
relation and the finite size scaling relation Eq. (14). Due to its better precision, we 
used above the lattice data obtained on lattices with four points in the temporal 
direction. It is expected [15] for quenched QCD that scaling may have set in for 
such lattices but the asymptotic scaling is violated by ~ 40-50 %, causing at most a 
factor of 2 change in Eq. (17). The systematic error in the finite size scaling relation 
of Eq. (15) due to this is more difficult to estimate but one does find similar results 
for larger temporal lattices than those used in ref. 19, suggesting again changes 
by factors of 0(1). Indeed, within the quenched QCD approximation it appears 
difficult to get away from the fact that the size of the universe at Tc measured 
in the units of the physical correlation length is essentially infinite, yielding an 
almost null result for the supercooling parameter. 

The real universe of course has dynamical fermions and one must therefore ask 
what happens in the case of full QCD. Unfortunately, in this case one has to worry 
about more input parameters in addition to the usual technical problem of putting 
fermions on the lattice. Assuming the current state of the art results, at best a 
very weak first order phase transition may be indicated for the realistic two light 
flavours and one heavier flavour case [20] . No signs of a flrst order transition were 
found on a lattice of linear size 4/Tc, suggesting a still larger correlation length 
at Tc, if flnite at all. A larger correlation length at Tc will increase the possibility 
of appreciable supercooling since a weaker phase transition will presumably yield 
broader region of metastability. Furthermore, a weaker phase transition will have 
lower latent heat and surface tension. These, in turn, may reduce 7 and therefore 
also the required degree of supercooling, -qh-, to achieve the desired mean free 
path. Thus, the scenario proposed in refs. [1-4] becomes more feasible for a very 
weak flrst order phase transition. However, it still looks rather improbable that 
the supercooling necessary for hadronic bubble creation will be achieved unless 7 
for the full QCD is lower by several orders of magnitude from its value for the 
quenched QCD used above. Whether future results from simulations of full QCD 
will still allow a small probability for homogeneous nucleation thus remains to 
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be seen. The alternative and more likely possibility is heterogeneous nucleation 
caused by 'impurities' left over from an earlier epoch. 
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